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In this letter we present a coherent picture for the evolution of Higgs mode in both neutral 
and charged s-wave fermion superfluids, as the strength of attractive interaction between fermions 
increases from the BCS to the BEG regime. In the case of neutral fermionic superfluid, such as 
ultracold fermions, the Higgs mode is pushed to higher energy while at the same time, gradually 
loses its spectral weight as interaction strength increases toward the BEG regime, because the 
system is further tuned away from Lorentz invariance. On the other hand, when damping is taken 
into account, Higgs mode is significantly broadened due to coupling to phase mode in the whole 
BEG-BGS crossover. In the charged case of electron superconductor, the Anderson-Higgs mechanism 
gaps out the phase mode and suppresses the coupling between the Higgs and the phase modes, and 
consequently, stabilizes the Higgs mode. 


The experimental search for Higgs boson in particle 
physics has made remarkable progresses UM- On the 
other hand, Higgs mode has also generated considerable 
interest in condensed matter and cold atom systems. 
Early in 1980s’, Raman scattering experiment has re¬ 
vealed an unexpected peak in a superconducting charge 
density wave compound NbSe 2 Isl], which was later at¬ 
tributed to the Higgs mode [J, . Signal of Higgs mode 

has also been observed in antiferromagnet TlCuCb by 
the neutron scattering , and recently in superconduct¬ 
ing NbN sample by terahertz pump probe spectroscopy 
in a nonadiabatic excitation regime [ 3 , . In cold atom 

system, Higgs mode has been observed near the super¬ 
fluid to Mott insulator phase transition of bosonic atoms 
in optical lattices at integer filling ili3- 


Theoretically, the simplest field theory where Higgs 
mode emerges is a relativistic U{1) field theory with 
Lorentz invariance in the symmetry broken phase. This 
occurs, for ex amp le, in the weak coupling BCS super¬ 
conductor ll|, [l2| or in the Mott-superfluid transition of 
Bose-Hubbard model at integer filling [HQ. However, 
in most condensed matter systems, Lorentz invariance 
only emerges with fine tuning and the generic symmetry 
is usually Galilean [1^. Thus, it is an interesting ques¬ 
tion to investigate how the Higgs mode evolves as the 
system is tuned away from the Lorentz invariance point. 
Moreover, in condensed matter systems, further compli¬ 
cations often occur because the Higgs mode is usually 
coupled to other elementary excitations which leads to 
its damping [HSi- In this Letter, we investigate these 
issues in the context of the BEC-BCS crossover model. 
In the BCS limit, the system obeys approximate Lorentz 
symmetry due to particle-hole symmetry and is expected 
to host Higgs mode. In the BEG limit, it is a condensate 
of molecular bosons and obeys the Galilean invariance. 
It thus provides a unique system to describe the fate of 
Higgs mode as the system is tuned away from Lorentz 
invariant limit. In addition, due to tunable interactions 





FIG. 1: (Golor online) v'Aq/u' and u!' jv! as functions of the 
scattering length (j — l/kpOs. In the inset we show v” Aq/u" 
as a function of C,. 


in the BEC-BCS crossover, it also provides a great plat¬ 
form to investigate the interaction effects on the Higgs 
mode due to coupling to collective and quasi-particle ex¬ 
citations [ 2 ^. 

We investigate these questions based on the time- 
dependent Ginzburg-Landau formulation of the BEC- 
BCS crossover, 

S = j dtSyi[(l)*{-iudt + vd'^ - (1) 


where (j) is the Ginzburg-Landau order parameter. The 
various parameters u,v,r,b and m* can be computed 
along BEC-BCS crossover in terms of the chemical poten¬ 
tial n, temperature T and C = l/(fcFas), where as is the 
s-wave scattering length. Within the Nozieres-Schmit- 
Rink 2l| framework, this can be calculated as detailed 
in the supplementary material 22|. The coefficients of 
the time derivative terms u = u' + iu" and v = v' + iv" 
are complex in general. The real parts u' and v' describe 
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the propagating behavior of the cooper pair field, while 
the imaginary parts u" and v" describe its damping due 
to coupling to the fermionic quasi-particles. A plot of 
various parameters are given in FigHJ We note the fol¬ 
lowing features. 

(i) Consider the real parts u' and v' in the BEC-BCS 
crossover. In the BCS limit, u'/ v'Aq —^ 0 because of 
the approximate particle-hole symmetry in the weak- 
coupling BCS theory while Aq = \/r/b is the mean field 
value of order parameter. As a result, the system ac¬ 
quires an emergent Lorentz invariance, and one expects 
the emergence of Higgs mode, together with the standard 
Anderson-Bogoliubov mode for neutral fermion super¬ 
fluid. In the BEC limit, however, v'Aq/u' ^ Aq/\^\ <C 1, 
and we can neglect the z;'-term. This leads to a Galilean 
invariant neutral boson theory, for which only Bogoliubov 
mode exists. 

(ii) The damping terms {u") becomes important as 
one moves to the BCS side, because of the decreasing 
fermionic excitation gap and as a result, a stronger cou¬ 
pling of the pairing field to the quasi-particle excitations. 
This corresponds to finite lifetime of Cooper pairs at fi¬ 
nite temperature. We will show that the damping u"- 
term itself will generate considerable effect for the ap¬ 
pearance of the Higgs mode different from that in a pure 
Lorentz invariance theory. In the BEC limit, the imagi¬ 
nary parts vanishes within NSR. On the other hand, we 
find that whenever they are nonzero, v"A q/u" 1 for 
the entire crossover regime and we shall thus neglect v"- 
term altogether in the following discussion. 

Spectral Weight Transfer without Damping. To inves¬ 
tigate the evolution of Higgs mode as the system is tuned 
gradually from its Lorentz-invariant BCS limit towards 
the Galilean invariant BEC limit, we shall first neglect 
the damping terms in EqUand study the transfer of spec¬ 
tral weight between the Higgs and Goldstone modes. In 
the symmetry broken state, we can write the order pa¬ 
rameter (j) = Aq -b da -b idp, where do and dp describe 
amplitude and phase fluctuations, respectively. In terms 
of da and dp and with u" = v" = 0, we can write the 
action Eq. |T]in the Fourier space as 


with <I>(w, k) = 
given by 



(da(w, k), dp(w, k)) and the kernel G is 


—iu'uj 


2r 


lU U) 


-b 


( 3 ) 


with k = |k| and ^k = k^/2m*. Two branches of spec¬ 
trum can be identified, with mode frequencies given by 



In the BCS limit, v'Aq ^ u' and solutions can be writ¬ 
ten as uj-{k) = k/y/2m*v' and a;+(fc) = a/(^ k + 2r)jv'- 



0 

0.5 1 

0 

0.5 1 

0 

0.5 1 




kf 


k^ 

« 10^ 

{a2) 

10^ 

(b2) 

1 

10^ 

(c2) 

1 


1 


1 


0 2 4 

0 2 4 

0 2 4 


oj/Aq 


oj/Aq 


Oj/Aq 


FIG. 2: (Color online) Spectral function Aao(k, ai) in the ab¬ 
sence of damping term. Aaa(k,ca) as a function of k (in unit 
of 1/5) and uj (in unit of Aq) for three different interaction 
strength 5 = —l/ikpas), C = —7 for (a), f = —3 for (b) 
and 5 = “1 for (c), corresponding to different gaps Aq/Ef = 
10“®, Aq/Ef = 4X 10“® and Aq/Ef = 7x 10~^, respectively. 
(a2-c2): Aaa(k,a;) as a function of u) for k = 0.1/5 (purple 
dashed line) and k = 0.01/5 (blue solid line). T/Tc = 0.9 and 
S is taken as 10“"^Aq. 


the first being the Goldstone mode with linear dis¬ 
persion and the second Higgs mode, with Higgs gap 
a;+(0) = \j2r/v' = \j2hA/f/v'. Using the facts that 
v' = 70^(/{‘i)vQ/ and b = 7pS in the BCS 
limit, one finds a;+(0) = 2Ao, as expected for a Lorentz- 
invariant theory. Here [3 = l/ksT is the inverse temper¬ 
ature, C,{n) is the Riemann-Zeta function and z/q is the 
density of state at the Fermi energy cp- 

In the BEC limit, u' ^ v'Aq, we find uj-{k) = 
\/Ck(Ck + 2,r)/u' is the Bogoliubov mode while the other 
mode 0 Jjr{k) = 02f]^/v' + 2r/v' -b (u'/v'Y has a gap 
~ |/i|, of order of binding energy of the molecule in the 
BEC limit. The existence of the gapped mode is a reflec¬ 
tion of the fact that our bosonic field ^ is a composite of 
two fermions and disappears in the infinite binding limit 
where only Bogoliubov mode exists as it should. 

In between these two limits, Lorentz invariance is bro¬ 
ken and the coupling between the amplitude and phase 
degrees of freedom becomes stronger, as characterized 
by the off-diagonal term iu'ui. We note that for low en¬ 
ergy Bogoliubov excitations, such coupling is small, but 
for gapped Higgs mode, it provides significant mixing of 
the amplitude and phase. To characterize such mixing, 
we calculate the spectral function for the amplitude 5a, 
given by Aaa{K^) = -:^IniC/aa(k,w-b id). Explicitly, 
this can be written as 

Aaa(k, w) = A+(k)d(a; - w+(k)) -b A_(k)d(w - w_(k)) 

( 5 ) 

with A+(k) and A_(k) being the spectral weight den- 
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sities associated with two modes lo- and a;+ |22l| . In 
FigH] (a,b,c), we plot the spectral function Aaa(k, w) for 
three representative values of C (corresponding to differ¬ 
ent Aq/Ep). Two features can be noticed immediately. 
First, the Higgs gap increases beyond 2Ao of the BCS 
limit as interaction strength increases. Secondly, there is 
increasing spectral weight transfer from the gapped Higgs 
mode to the gapless mode. One can show explicitly that 
A+/A- = Av'r‘^{u''^ /2m*{k‘^/2m* + 2r))~^, which 
indicates the gradual increasing of the mixing between 
phase and amplitude degrees of freedom. 

Including Damping Term. Due to the presence of 
damping term, the time-dependent Ginzburg-Landau 
theory is not a pure Lorentz invariant C/(l) theory. Thus, 
at any finite temperature, even in the BCS limit, the peak 
of Higgs excitation uj+ (k) will not as sharp as discussed 
above. To calculate the equilibrium spectral weight in 
the presence of damping, we need to introduce the so- 
called Langevin force which satisfies the follow¬ 

ing conditions, (r]{t',x.')r]{t,x.)) = (? 7 *(t', x') 7 y*(t, x)) = 0 , 
and {ri*(t',x.')r]{t,E)) = 2u''kBT6(t — t')d{:x. — :x.'). In¬ 
cluding the corresponding term in the action as S'l = 
J dtd^x{(j)*r] + (/rf), we obtain the equations of mo¬ 
tion for 5a and 5p, by setting d{S + SL)/dSa = 0 and 
a(5 + 5L)/a<5p = 0 , 

[—vuP -t- Ck + 2r) 5a — iuujSp + 7]' = 0, (6) 

+ ^k) 5p + iuuj5a + rj" = 0, (7) 

where p' and rj” are the real and imaginary parts of the 
Langevin force p, respectively. The spectral functions for 
the amplitude fluctuation is given by, using fluctuation 
dissipation theorem, 

_u"uj I — -I-Ckp + 

“ “ 2 I - (uw)2 + (-ua;2 + + ^k + 2r)|2 ' 

( 8 ) 

By comparing Fig. [3] with Fig. [3J one can see three 
important features brought about by including the damp¬ 
ing term. First, the spectral weight transfer is enhanced. 
For instance, for = —7, there is almost no spectral 
weight transfer in the absence of damping (Fig. UKa)) 
while in the presence of damping, for very small k <C 1 /^, 
Haa(k, w) exhibits a clear peak at the energy of Bogoli- 
ubov mode, with a weight proportional to u” [2^. Sim¬ 
ilar enhancement of spectral weight transfer can also be 
easily seen in Fig. [3Kb) for (/ = —3. Secondly, also for 
A: <C 1/^, in the BCS limit, the location of Higgs peak is 
substantially reduced from ^2r/v' to ^J2r/v' — u"'^/v'"^, 
as shown for (/ = —7 and —3 in Fig. jSKa) and (b), re¬ 
spectively [ 2 ^. Thirdly, as k starts to derivate from zero, 
the Higgs mode quickly loses its identity, due to strong 
hybridization with the Bogoliubov mode. For instance, 
even for fc = 0 . 1 /^, as displayed by the purple dashed line 
in Fig. |3Ka2-c2), no feature of sharp peak is observed in 
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FIG. 3: (Color online) Spectral function Taa(k,a;) in presence 
of damping term. Aaa(k,a;) as a function of k (in unit of 1/^) 
and LJ (in unit of Aq) for three different interaction strength 
C = —l/(kFas), C = —7 for (a), ^ = —3 for (b) and ^ = —1 for 
(c), corresponding to different Aq/Ef = 10“®, Aq/Ef = 4 x 
10 “® and Aq/Ef = 7x 10“^, respectively. (a2-c2): Haa(k,a;) 
as a function of oj for k = 0.1/^ (purple dashed line) and 
k = 0.01/if (blue solid line). T/Tc = 0.2 and 5 is taken as 
10"^ Ao. 


-^aa (k,w). And for C, = —1, no sharp peak exists even 
for k = 0.01/^. 

Effects of Coupling to External Gauge Fields. Now we 
understand that, in the weakly interacting BCS side of a 
neutral superfluid, the appearance of Higgs mode suffers 
significant broadening due to finite M"-term at finite tem¬ 
perature, which couples the Higgs mode to the collective 
Bogoliubov excitations. Therefore, if we further consider 
the presence of coupling to external electromagnetic field 
for the case of charged fermions, the Bogoliubov mode 
is gapped out by the Anderson-Higgs mechanism. Thus, 
we expect that the Higgs mode is easier to observe in 
the charged case. To incorporate this effect of external 
electromagnetic field, we introduce the gauge potential 
ip(t, x) and extend the action as 


5-0 = 


J dtd^x{(l)*[—iu{dt — 2eip) -\- v{dt — 2eipY 


(9) 


b. 

- r]d> -\— 

2 m* ' 2 ' 




where e is the charge of the electron. Following the same 
procedure as before, we find that the coupling between 
5a and 5p is modified and is now proportional to k^ 


2uu)k'^ 

-b 327rue^AQ 


(5a(w, k)5p(-a;, -k). 


( 10 ) 


As a result, the original gap less phase mode is gap ped to 
a finite frequency, w(k) = + 167re^AQ/m*, which 

is known as the Anderson-Higgs mechanism. Thus, the 
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FIG. 4: (Color online) Spectral function j4aa(k, oj) for the 
charged case. Aaa(k, cu) as a function of k (in unit of 1/^) 
and OJ (in unit of Ao) for three different interaction strength 
C = —l/ikpas), C = —7 for (a), ^ = —3 for (b) and ^ —1 for 
(c), corresponding to different Ao/E f = 10“®, Aq/Ef = 4 x 
10“® and Aq/Ef = 7x 10“^, respectively. (a2-c2): Aaa(k,a;) 
as a function of cn for k = 0.1/^ (purple dashed line) and 
k = 0.01/if (blue solid line). T/T^ = 0.9 and 5 is taken as 
10"^ Ao. 


large energy separation between this gapped phase mode 
and Higgs mode strongly suppresses their coupling. A 
further consequence of the modification is that at long 
wave length fc —>■ 0, the coupling between phase and 
amplitude mode becomes small. The spectral function 
Aaa(k, w) for charged case is plotted in Fig. 0] In 
sharp contrast to neutral case Fig. 01 the presence of 
damping term has almost no effect on Higgs mode, and 
there is always a peak located at w = 2Ao. In this case, 
as attractive interaction increases and the system gradu¬ 
ally loses its Lorentz invariance, the peak becomes more 
and more broad. 

Conclusion. In summary, we have investigated the evo¬ 
lution of Higgs mode in the BEC-BCS crossover for both 
neutral and charged Fermi superfluid. Our main con¬ 
clusions include: i) Towards the BEG side, as the sys¬ 
tem gradually loses the Lorentz invariance, the Higgs 
mode is pushed to very high energy and the spectral 
weight is transferred to Bogoliubov mode, ii) In the BCS 
side, damping terms arises in the Ginzburg-Landau the¬ 
ory, due to coupling between Gooper pair field and the 
fermionic quasi-particles, and strongly couples the Higgs 
mode to the gapless phase mode in the neutral superfluid, 
which enhances the spectral weight transfer and washes 
out features of Higgs mode at finite momentum, (iii) 
For the charged case, the phase mode is gapped out by 
coupling to external electromagnetic field, and the Higgs 
mode becomes much more stable. 

Our results also deepen our understandings of Higgs 
mode in superconductor. The physical picture behind 


the observation of Higgs mode in a BGS superconduc¬ 
tor is much more subtle and its observability is not 
merely guaranteed by Lorentz symmetry. While the 
damping terms broadens the Higgs peak, the Anderson- 
Higgs mechanism alleviate the coupling between Higgs 
and phase mode and as a result, Higgs mode remains 
at energy 2Ao. As for cold atom system, because of the 
cooling limit, so far we can not reach Fermi superfluid for 
C < — 1. However, our results show no Higgs feature in 
spectral function for C > — 1. On the other hand, with re¬ 
cent development of synthetic gauge field, there are many 
proposals to generate a synthetic dynamic gauge field in 
cold atom system [2^ . If such a dynamic gauge field can 
be experimentally realized and coupled to fermions, the 
Anderson-Higgs mechanism will be activated and a Higgs 
mode will be observed. This can be used as a way to test 
our theory. 
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I 

SUPPLEMENTARY MATERIALS 


Time-dependent Ginzburg-Landau theory of BEC-BCS crossover 


A time-dependent Ginzburg-Landau theory can be constructed for the entire BEC-BCS crossover in the vicinity of 
Tc [II- The partition function takes the form 2 = J with 




( 1 ) 


where tpa- are Grassman fields and g is the contact interaction between fermions of opposite spins, g is the chemical 
potential which is determined by requiring the number density to be equal to n. To investigate the fluctuation effects 
in the Cooper channel, we use a Hubbard-Stratonovich transformation to decouple the interaction term in the Cooper 
channel and then integrating out the fermions. We obtain an effective theory for the bosonic field A(t, x), which 
represents the cooper pair field. Straightforward calculations yield the partition function in terms of field A as 
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i4(A, A) exp 
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is the Gor’kov Green function. 

In the vicinity of the phase transition the gap parameter A is small and an expansion in terms of A becomes 
possible. Including both the spatial and time derivatives (after Wick rotation) and retaining the parameter A up to 
the forth order we obtain an effective action as 
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where u = u' iu" and v = v' iv" are complex in general and all the parameters can be expressed in terms of 
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microscopic parameters as 
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( 11 ) 


In the above equations, -/V(a) = l/(exp(/3a) +1) is the Fermi distribution function and a = ^k —/r with Ck = k'^/2m. 
Function 0(2/x) is the heaviside step function. Explicitly, the parameter b is the result of one-loop calculation with 
four fermion propagators 


1 f d^k 1 1 

/32 J ( 27 r)^ {—iijJn + fc2/2m — /x)2 (iojn + A:2/2m — /i)2 


( 12 ) 


The other parameters u, v, and r are all derived from the inverse vertex function F ^(a;„,k), which after the 
standard renormalization by replacing g with the two-body scattering length Og, is given by 


F ^(w„,k) = 


m 

dTTOs 


f d^k 

J Wr 


1 - jV(ek - m) - -^(Ck-q - +) 

— iuJn + Ck + Ck-q ~ 2/r 



(13) 


To derive the time-dependent Ginzburg-Landau equation, we first analytically continue vertex function to real fre¬ 
quency iuJn —>■ w + 10+. This procedure generates a time-dependent term with parameter u and v. The detailed 
derivation is as following. 

The frequency dependent part of r“^(w,k) is 
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Then we expand it in series of small oj as 
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We define the parameters as u = f '^ = f ( 2+3 both can be calculated by 

contour integration. 
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( 16 ) 


where c denotes the contour in Fig. [T] There are infinite first-order poles Zn = fJ- + and one second order 

pole Zjj = ^ + rji. The contour integration can be evaluated in the summation of the residuals as 



FIG. 1: The contour “c” in the calculation of Eq. (1161) . The dots denote the first-order poles Zn = and the 

cross “x” denotes the second order pole Zn = fi + rji. 
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7r2 2^ 
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n— — C)0 
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(17) 


Calculation shows that J2n=-oc pure imaginary due to the symmetry of the z„ pole 

locations with respect to the horizontal axes. Then it can be written as X)n2_oo '^^'*"( 2 ”+!)^^^''^^ “ 


v2*e::o 


-Hoo A/Vl-K(2n-H)7r//3;r)^-Sgn()() 


(2n+iy 


. Hence, the parameter u is calculated as 


(2m)^/^ 

IGtt^ 


'2 V2Py/\iA - ~ 

TT < ^ 


n—0 


{2n+iy 




(18) 


In the same manner, the parameter v can also be calculated as shown in Eq. (6) and (7). We should note that while 
the expressions for u and others look different from the standard expression, as given in ref. [l|, they in fact reduce 
to the same expressions. We found that this form is more convenient to use the above expression when dealing with 
higher order time-derivative terms. 

In the BCS and BEC limits all the parameters can be analytically derived as shown in Table I. 


Spectral weight function in the case without damping term 


If we ignore the damping term by taking u” = 0 the action can be written as 

duj - 


S = 


-$(-a;,-k)a-^$(w,k) 


27r (27r)3 


(19) 
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Parameters 

BCS limit 

BEC limit 

u' 

0 

7rn{ep) 

u" 

rb 

00 

0 

v' 


7ra(ep) 


64yeF|Mp/2 

v" 


0 

1 

2m* 

1 7i/(eB)eF Atn\ 

2m 12-K2(kBT)2 Sloj 

1 7riy(€^) 

r 

iz(eF) in ^ 

2 y2,/rF \/2|hl) 

b 

/-/o^ 

87r^(fcgT)^ / 

TTV{ep) 

32,/eF|Mp/2 


TABLE I: Asymptotic behaviors of the parameters in the time-dependent Ginzburg-Landan theory in the BCS and BEG limits. 


with $(u;, k) = ((5a(w, k), Sp(u!, k)) and the kernel Q is given by 


= 


—v'uP' -b + 2r m'w 

—iu'uj + Ck y ’ 


with k = |k| and ^k = /2m*. Then the amplitude-amplitude correlation function can be easily calculated as 


Qaaiy^ 1 h) — 


+ ^k 


+ {—v'uj^ + ^k)(—+ ?k + 2r)' 
Straight forward calculation yields the spectral function as 

Aaa(w,k) = --ImGaaiu} + iS,k) 

= A+(k)5(u; — cu+(k)) + A_(k)(5(cu — u;_(k)), 

where the mode frequencies are given as 


and the spectra weight density 


- ~/;r- + ^^\l772+TJ4+ 773 + r) 


A,ik)= 


A_(k) = 


2v'^uj+{col - L 0 ^_)’ 
-v'uj'i + ^k 


2v'^uj-{loI - UJ^_)' 

At BCS limit the ratio of the two spectral weight densities can be approximately calculated as 


A_(k) y^k'^/2m*{k‘^/2m* -I- 2r) 

A+(k) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


At BCS limit we have u'/v' —>■ 0, this ratio vanishes. This spectral weight transfer is shown in Fig. 2 in the main 
text. 


Spectral weight function in the case with damping term 

The spectral weight function of the amplitude mode in the case with damping term u" is 

Ann. 
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where the eigen mode frequencies are 

+ r u''^ — 
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2?;'2 
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„/3 
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For small momentum they can be approximated as 


2rek 


2v'r + m' 2 — u"2 ’ 


w+ = 


12v'r + u'2 - m "2 2v'r + 2m' 2 - 2u"2 


?k- 


(28) 


vf 2 v'(2v'r + u'2 — lt"2) 

Compared with the case without damping term we see that the gap of the Higgs mode is reduced from ^y2rjv' to 
yj2rlv' — m"2/u' 2 at BCS limit. 

For small ^k the spectral weight on the Goldstone mode can be calculated as 


Haa(w_,k) = 


8w'2w_ ■ 


(29) 


Different from the case without damping term, we see that in the case with damping term the spectral function has 
a weight proportional to v!' on the Goldstone mode. 


The spectral weight function in the case with Coulomb interaction 

A time-dependent Ginzburg-Landau theory with Coulomb interaction can be cast as 0 


F = 


[ dt(fxl --I- Af — iu{dt — 2e(/)) -I- v(dt — 2e()>)^-^-r) A -|- -AAAaV 

7 1 Stt V 2m* J 2 J 


(30) 


where e is the electric charge and x) is the electric field. By taking a symmetry breaking A —>■ Aq -f i5a -I- i6p we 
can have a free energy for the low energy excitations in the momentum space as 

/ du] f ^2 

^ (271^ 1 k)(5p(-w, -k) -I- da(-w, -k)(-uw^ -I- -I- 2r)SaiuJ, k) -|- Sp{-uj, -k) 


(—-I- -— -)Sp{u}, k) -I- 4meAo(/)(—w, — k)da(w, k) — 4veujAo4>{—oJ, —k)(5p(w, k) -|- due^Ap^)— w, —k)(/)(a;, k) j. 

(31) 


We integrate out the electric field </) and obtain 
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Then the spectral functions can be calculated as 


(32) 
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(33) 
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